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Abstract. We present J =2-1, J=3-2, J=4-3 12 CO maps as well as 3=1-1, J=3-2 13 CO and 492 GHz [CI] 
measurements of the central region in M 51. The distribution of CO is strongly concentrated towards the spiral 
arms. The center itself is poor in, though not devoid of, CO emission. The observed line intensities require 
modelling with a multi-component molecular gas. A dense component must be present (n( Ha) ~ 10 3 ) with kinetic 



temperature Tkin ~ 100 K, combined with either a less dense (f» 10 cm" ) component of the same temperature, 
or a more dense (n(H2) « 3 x 10 3 cm -3 ) and much cooler (7km = 10-30 K) component. Atomic carbon amounts 
are between 5 and 10 times those of CO. Much of the molecular gas mass is associated with the hot PDR phase. 
The center of M 51 has a face-on gas mass density of about 40 ± 20Mq pc~ 2 , and a well-established CO-to-H2 
conversion ratio X four to five times lower than the standard Galactic value. 
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1. Introduction 

Molecular gas is a major constituent of the interstellar 
medium in galaxies. This is particularly true for star- 
forming complexes in the spiral arms, but strong cir- 
cumnuclear concentrations of molecular gas are frequently 
found also in the inner kiloparsec of spiral galaxies. We 
have observed a sample of nearby spiral galaxies in various 
CO transitions and in the 492 GHz 3 Pi- 3 P [CI] transition 
in order to determine the physical condition of molecular 
gas in their inner parts. Results for NGC 253 (Israel et 
al. 1995), NGC 7331 (Israel & Baas 1999), NGC 6946, 
M 83 = NGC 5236 (Israel & Baas 2001) as well as IC 342 
and Maffei 2 (Israel & Baas 2003) have been published. In 
this paper, we present results obtained for the interacting, 
two-armed spiral M 51 = NGC 5194 (see TableHJ). 

M 51 was one of the first galaxies mapped in CO 
line emission, and it has been observed many times since 
(Rickard et al. 1977; Scoville & Young 1983; Rydbeck et 
al. 1985; Sandqvist et al. 1989; Garcia-Burillo et al. 1993a, 
b; Berkhuijsen et al. 1993; Nakai et al. 1994; Young et al. 
1995; Mauersberger et al. 1999; Wielebinski et al. 1999; 
Dumke et al. 2001; Paglione et al. 2002; Tosaki et al. 
2002), including mapping observations with millimeter ar- 
rays (Lo et al. 1987; Rand & Kulkarni 1990; Rand 1993; 
Scoville et al. 1998; Aalto et al. 1999; Sakamoto et al. 1999; 
Regan et al. 2001). 



Table 1. 



Galaxy parameters 



M 51 



Type" 

Radio Center : 
R.A. (B1950) 6 
DecL^^O) 6 
R.A. (J2000) 6 
Decl.(J2000) 6 
V c ' 

LSR 

Inclination i c 
Position angle P c 
Distance D d 
Scale 



SA(s)bcp 

13 h 27 m 46.3 s 

+47°27'10" 

13 ,l 29 m 52.7 s 

+47°11'42" 

+464 kms" 1 

20° 

170° 

9.7 Mpc 

21.3 , '/kpc 



Notes to Table 1: a RSA (Sandage & Tammann 1987); b Turner 
& Ho (1994); c Tully (1974) d Sandage & Tamann (1975) 



2. Observations 

All observations described in this paper were carried out 
with the 15m James Clerk Maxwell Telescope (JCMT) 
on Mauna Kea (Hawaii) 1 . At the epoch of the mapping 
observations (1990-1996) the absolute pointing of the 
telescope was good to about 3" r.m.s. as provided by 
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Table 2. 
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Observations Log 



Transition 


Date 


Freq 




Beam 




t(int) 




Map Parameters 










Size 






Points Size 


Spacing 




(MM/ Y Y) 


(GHz) 


(K) 


( ) 




(sec) 




(") 


\ ) 


12 CO J =2-1 


12/90; 05/91 


230 


920 


21 


0.69 


480 


40 


30 x 50 


10 


CO J =3-2 


12/93 


345 


1475 


14 


0.53 


400 


18 


40x50 


12 




05/97- 05/01 




660 




0.60 


70 


400 


60x120 


6 


12 CO J=4-3 


03/94; 04/96 


461 


3950 


11 


0.50 


600 


16 


24x20 


6 




12/01 




2100 




0.52 


1200 


6 


18x12 


6 


13 CO J=2-l 


06-95 


220 


370 


21 


0.69 


2520 


1 






13 CO J=3-2 


05-97; 05/01 


330 


600 


14 


0.63 


3600 


9 


15x15 


7.5 


CI 3 Pi- 3 P 


11-94 


492 


6675 


10 


0.53 


960 


4 







pointing observations with the JCMT submillimetre 
bolometer. The spectra were calibrated in units of 
antenna temperature TjJ, correcting for sideband gains, 
atmospheric emission in both sidebands and telescope 
efficiency. Calibration was regularly checked by obser- 
vation of a standard line source. Further observational 
details are given in TableEl Most of the observations were 
carried out with the now defunct receivers A2, B3i and 
C2. Observations in 2001 were obtained with the current 
receivers B3 (330/345 GHz) and W/C (461 GHz). Full de- 
tails on these receivers can be found at the JCMT website 
( IhttpV/docs.jach.hawaii.edu/JCMT/HET/GUIDETI l. 
Note that the angular beamsizes used correspond to 
linear resolutions ranging from 470 to 1000 pc at the 
distance assumed for M 51. Up to 1993, we used a 
2048 channel AOS backend covering a band of 500 MHz 
(650 km s" 1 at 230 GHz). After that year, the DAS 
digital autocorrelator system was used in bands of 500 
and 750 MHz. Integration times (on+off) given in Table|3 
are typical values appropriate to the maps. Because 
of the relatively good weather conditions and its very 
close sampling, the 345 GHz J=3-2 12 CO map (shown 
in Fig.EJ) should be considered the most reliable. Both 
the 230 GHz J=2-l 12 CO and the 461 GHz J=4-3 
12 CO maps were obtained under less favourable weather 
conditions (as evidenced by the system temperatures 
listed in Table0). Features in the 230 GHz map show 
offsets from their counterparts in the 345 GHz map by as 
much as 8", and the peak in the 461 GHz map is offset 
from its 345 GHz counterpart by 5.5". As these map 
position differences are less than three resp. two times 
the r.m.s. pointing accuracy, we do not consider them to 
be physically significant. 

M 51 is observed almost face-on, lines are narrow and 
sufficient free baseline was available to subtracted second 
or even third order baselines from the profiles. Finally, all 
spectra were scaled to a main-beam brightness tempera- 
ture, T m b = Tj^/r] m b; values for rj^ are given in TableEl 

Much of the following analysis is based on line ratios 
at two specific positions A and B. Position A is that of 
the (radio) nucleus, position B is centered on a large and 



bright CO complex (called Giant Molecular Association 
1 by Tosaki et al. 1994) in one of the inner spiral arms 
at a radial distance of 0.9 kpc northwest of the nucleus. 
Measurements of these positions are much more accurate 
than those of individual map positions. They were usually 
observed more than once, with better pointing (typically 
1.5" — 2" r.m.s.) and with significantly longer integration 
times resulting not only in higher signal-to-noise ratios 
but also in better baselines and higher position reliabil- 
ity. Relatively accurate line intensities for full-resolution 
profiles at the two positions are summarized in Tabled 
In this Table, we also give extrapolated line intensities at 
other resolutions, derived as follows. We compared inten- 
sities at the two positions in maps at full-resolution and 
in maps convolved to the desired resolution. The accurate 
full-resolution measurements in Tabled were then multi- 
plied by the empirical scaling factors thus determined. We 
verified that this is a robust procedure. Pointing errors do 
result in different intensities at any given position, but 
scaling factors change much more slowly. 

3. Results 

3.1. CO distribution 

The distribution of molecular gas as traced by CO shows 
a minimum at the position of the galaxy nucleus (position 
0", 0" in Fig. CO is strongly concentrated towards the 
major spiral arms. The CO maps in Fig.[5] clearly show 
the outline of the inner parts of the western spiral arm, 
as well as part of the inner eastern arm. In all observed 
transitions, the emission peaks coincide with the brightest 
CO peaks seen in the J=l-0 maps (see e.g. Scoville et al; 
1998 or Regan et al. 2001) which are associated with star- 
forming regions. 

The north-south position- velocity maps (Fig.EJ), par- 
alleling the major axis of M 51, show emission covering a 
velocity range of 150 kms^ 1 , which becomes 440 km s -1 
after correction for inclination. The velocity distributions 
east and west of the major axis are practically identical, 
at least at our resolution. Within a radius R < 30" the 
molecular gas is in rapid solid-body rotation. Again, the 
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R.A. offset(arcsec) from 13 27 46.28 



R.A. offset(orcsec) from 13 27 46-31 



R.A. offset(arcsec) from 13 27 45.31 



R.A. offset(orcsec) from 13 27 46.29 



Fig. 1. Contour maps of emission from M 51 integrated over the velocity range Vlsr = 300 to 600 kms 1 . North is 
at top. Offsets are marked in arcsec relative to the M 51 nucleus at (0,0) with B1950 coordinates given corresponding 
to J2000 coordinates a = 13 h 29 m 52. s 7, 5 = +47°ll'43" (cf. Tabled . Left to right: CO J=2-l at 21" (985 pc) 
resolution, CO J=3-2 convolved to 21", CO J=3-2 at 14" (656 pc) resolution and CO J=4-3 convolved to 14" 
resolution. Contour values are linear in f T m bdV. Contour steps are 4 Kkms -1 (2-1 full resolution, 3-2 convolved, 
4-3 convolved) and 8 Kkms -1 (3-2 full resolution) and start at 0. 



Table 3. 



Central CO and CI line intensities in M 51 



Transition 


Resolution 


Tmb / T mb dV 


T mb / T mh dV 




(") 


(mK) (Kkms- 1 ) 


(mK) (Kkms -1 ) 






Position A: 


Position B: 






Center 


NW Arm GMA-1 



12 CO 


J=2-l 


21 


260 


38.8±8.5 


805 


54.0±7.5 






43 




44.5±9.0 




36.0±7.2 


13 co 




21 






77 


7.5±1.0 


12 co 


J=3-2 


14 


375 


42.7±6.5 


1535 


64.5±9.5 






21 




30.2±6.0 




55.0±9.5 


13 co 




14 


28 


3.3±0.7 


43 


5.5±1.5 


12 co 


J=4-3 


11 


160 


27.0±4.5 


520 


52.5±7.5 






14 




23.0±4.5 




33.0±7.0 






21 




19.0±3.8 




24.0±4.0 


[CI] 


3 Pi- 3 Po 


10 






565 


28±5 






21 








(13±4) 



Note to Table 3: Position A: beam centered on nucleus; Position B: beam centered on Giant Molecular Association 1 in NW 
spiral arm at offsets Aa = —10", A8 = +15" with respect to the nucleus. 



relative lack of molecular gas in the center (position +464 
km/s, 0") is obvious, as is the strong molecular emission 
near the radius where differential rotation takes over. As 
the pV-maps of M 51 are virtually indistinguishable from 
the major axis pV-diagrams of other, often more highly 
tilted galaxies, they illustrate an interesting point. In those 
galaxies, the existence of a compact central feature with 
a steep linear velocity gradient is often interpreted as the 
signature of a circular-symmetric physical structure such 
as e.g. a rotating circumnuclear disk or torus. However, 
the almost face-on image of M 51 clearly shows the ab- 
sence of such large-scale structure and the diagrams in 
Fig-El simply reflect the strong molecular enhancement of 
the inner spiral arms. 



3.2. Line ratios 

Using our own and literature data, we have determined 
line intensity ratios of the observed transitions at both 
the center (position A), and at the position of the molec- 
ular cloud complex in the northwestern arm (position B). 
All 12 CO line ratios listed in TableH are for a 21" beam 
(roughly corresponding to a circle of 0.5 kpc radius). The 
isotopic intensity ratios 12 CO/ 13 CO also refer to that 
beam size in the J =1—0 and J=2— 1 transitions, but to 
a smaller 14" beam in the J=3-2 transition. However, the 
isotopic intensity ratios appear to vary little in the in- 
ner part of M 51 (cf. lit. cited) so that we expect that 
observation with a larger beam would have yielded very 
similar values for the J=3-2 transition. Note that the 
[CII]/CO(2-l) ratio was determined after convolution of 
the J=2-l 12 CO map to the 55" beam of the [CII] obser- 
vations. The [CII] intensities (Crawford et al. 1985; Stacey 
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Fig. 2. Position- velocity maps of CO J—i-2 emission from M 51 in position angle 10°. Left: east of the nucleus, 
averaged Aa = 20"-0"; right: west of the nucleus, averaged over Aa = 0"- -20". Contour values are linear in steps of 
T m b = 4 K, starting at step 1. Horizontal scale is Vlsr 

Table 4. Integrated line ratios in the centre of M 51 



Transitions Center NW Arm GMA-1 





Observed 


Modelled" 


Observed 


Modelled 


12 CO (l-0)/(2-l) 6 


1.25±0.10 


1.28 


1.40±0.20 


1.35 


12 CO (3-2)/(2-l) c 


0.78±0.23 


0.75 


1.02±0.16 


0.87 


12 CO (4-3)/(2-l) c 


0.49±0.10 


0.41 


0.44±0.08 


0.46 


12 CO/ 13 CO (l-0) d 


7±1.5 


6.7 


8±1.5 


8.2 


12 CO/ 13 CO (2-l) c ' d 


6±1 


6.9 


8.7±1.2 


8.1 


12 CO/ 13 CO (3-2) e 


12.9±1.8 


12.2 


12.5±1.7 


13.2 


CI/CO(2-l) c 


0.22±0.06 




0.24±0.08 




CII/CO(2-l) / 


0.46±0.09 









Notes: a. See Tabled] b. Ratios estimated from data by Nakai et al. 1994, NRO at 16" resolution; Garcia-Burillo et al. 1993, 
IRAM at 22". c. This Paper, JCMT at 21" resolution, d. Ratios estimated from data presented by Rickard & Blitz 1985, 
NRAO 65", Garcia-Burillo et al. 1993, IRAM at 22", Matsushita et al. 1997, NMA at 6"; and Paglione et al. 2001, FRAO at 
45" resolution, e. This Paper; JCMT at 14" resolution, f. From Crawford et al. (1985) and Stacey et al. (1991), KAO at 55" 
resolution; consistent with ISO flux (Negishi et al., 2001); JCMT convolved to 55" resolution. 



et al. 1991; see also Negishi et al. 2001) were converted to 
velocity-integrated temperatures to obtain the line ratios 
in TableU 

The 12 CO and 13 CO ratios of M 51 are quite different 
from those of the other galaxies (NGC 6946, M 83, IC 342, 
Maffei 2; Israel & Baas 2001, 2003) we have studied so far. 
The ratio of [CI] to J=2-l 12 CO intensities is relatively 
high, but that of [CII] to J=2-l 12 CO is not dissimilar 
from those seen in other galaxies. It resembles, in particu- 
lar, the ratios seen in galaxies with central starbursts. We 
suspect that this similarity is caused by the large (radius 
R Ki 1.3 kpc) aperture used to measure the [CII] line in 
M 51. We suspect that the very center of M 51, weak in 
CO and HI, contributes relatively little, and that most of 
the emission in the aperture comes from the inner spiral 
arms with their strong emission from star- forming regions. 
Bright [CII] emission indicates the occurrence of both high 
temperatures and high gas densities as the critical values 
for this transition are T^ n > 91 K and n > 3500 cm -3 . 



However, such values must be reconciled with the rela- 
tively high CO opacities implied by the modest isotopic 
intensity ratios in the lower CO transitions. In the arm po- 
sition B, isotopic ratios increase with increasing J level, 
but at the center position A the isotopic ratio is lowest in 
the J=2-l transition. This reflects significant differences 
in the molecular gas properties at the two positions sam- 
pled. 

4. Analysis 

4.1. Modelling of CO 

We have modelled the observed 12 CO and 13 CO line in- 
tensities and ratios with the large- velocity gradient (LVG) 
radiative transfer models described by Jansen (1995) and 
Jansen et al. (1994). These provide model line intensities 
as a function of three input parameters per molecular gas 
component: gas kinetic temperature Tk, molecular hydro- 
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Table 5. Model parameters 



Model 




Component 1 




Component 2 


Relative 




Kinetic 


Gas CO Column 


Kinetic 


Gas CO Column 


J=2-l 12 CO 




Temperature 


Density Density 


Temperature 


Density Density 


Emission 




T k 


n(H 2 ) iV(CO)/dV 
(cm" 3 ) (cm-^kms- 1 )- 1 ) 


Tk 


n(H 2 ) iV(CO)/dV 
(cm- 3 ) (cm-^kms- 1 )- 1 ) 


Component 




(K) 


(K) 


1:2 



1 (Center) 100 100 l.OxlO 17 60 1000 3.0xl0 17 0.65:0.35 

2 (Center) 150 1000 0.6xl0 17 20 3000 3.0xl0 17 0.55:0.45 

3 (Arm) 100 100 l.OxlO 17 150 1000 3.0xl0 17 0.80:0.20 



Table 6. Beam-averaged results 



Model 


Beam- Averaged Column Densities Face-on Mass Density Relative Mass 
iV(CO) N(C) JV(H 2 ) <r(H 2 ) a| as Components 1:2 
(10 18 cm- 2 ) (10 21 cm- 2 ) (M /p C - 2 ) 


Ratio A r (H 2 )//(CO) 
X 

(10 20 molcm- 2 /Kkms- 1 ) 


M 51; N K /Nc = 1250; ATHI) C = 2 x 10 20 cm" 2 


1 

2 
3 


0.36 4.4 
0.36 1.2 
0.42 3.3 


2.9 43 60 0.6 : 0.4 
1.0 15 20 0.4 : 0.6 
2.2 33 47 0.8 : 0.2 


0.74 
0.25 
0.41 



Note: a. Sum of H 2 and HI multiplied by 1.35 to take into account the contribution by Helium; b. See text, Sect. 4.2; c. HI 
column density valid for both Center and (nearby) Arm; Tilanus & Allen (1991). 



gen density n(H2) and the CO column density per unit 
velocity N(CO)/dV. By comparing model to observed line 
ratios, we may identify the physical parameters best de- 
scribing the actual conditions at the observed positions. 
Beam-averaged properties are determined by comparing 
observed and model intensities. In principle, with seven 
measured line intensities of two isotopes, properties of a 
single gas component are overdetermined as only five in- 
dependent observables are required. We found that no fit 
based on a single gas component is capable of matching 
the data at either of the observed positions. 

However, we obtained good fits based on two gas com- 
ponents. In order to reduce the number of free parameters, 
we assumed identical CO isotopical abundances for both 
gas components. In a small number of starburst galaxy 
centers (NGC 253, NGC 4945, M 82, IC 342, He 2-10), 
values of 40 ± 10 have been suggested for the isotopical 
abundance [ 12 CO]/[ 13 CO] (Mauersberger & Henkel 1993; 
Henkel et al. 1993, 1994, 1998; Bayet et al. 2004), some- 
what higher than the characteristic value of 20-25 for 
the Milky Way nuclear region (Wilson & Rood 1994). 
Although M 51 does not have a starburst nucleus but 
rather a low-luminosity AGN, vigorously starforming spi- 
ral arms already occur at small radii. For this reason we 
have adopted an abundance value of [ 12 CO]/[ 13 CO] = 40 
in our models. We identified acceptable fits by searching 
a grid of model parameter combinations (10 K < Tk < 
150 K, 10 2 cm- 3 < n(H 2 ) < 10 5 cm- 3 , 6 x 10 15 cm" 2 < 
N(CO)/dV < 3 x 10 18 cm- 2 ) for model line ratios match- 
ing the observed set, with the relative contribution of the 
two components as a free parameter. Solutions obtained 



in this way are not unique, but rather delineate a range 
of values in distinct regions of parameter space. For in- 
stance, variations in input parameters may to some ex- 
tent compensate one another, producing identical line ra- 
tios for somewhat different combinations of input param- 
eters. Among all possible solution sets, we have in any 
case rejected those in which the denser gas component is 
also hotter than the more tenuous component, because we 
consider the large pressure imbalances implied by such so- 
lutions physically implausible, certainly on the kiloparsec 
scales observed. 

The results for the two positions in M 51 are summa- 
rized in Table|3J The physical conditions applying to the 
molecular complex in the arm (position B) are very well 
determined by the observations. The emission is wholly 
dominated by fairly hot gas with a kinetic temperature 
of the order of Tki n = 100 K. Low-density gas with 
hh 2 " 100 cm- 3 contributes most (80 ± 10%) of the ob- 
served emission, the remainder (20 ± 10%) coming from 
gas at a moderately high density hh 2 " 1000cm- 3 (for- 
tuitously, the masses of the two components are in almost 
the same proportion - cf. TableEJ. These parameters pro- 
vide a very good fit as is evident from Table^J no other 
set of values does. 

In contrast, the physical conditions determining the 
emission from the center of M 51 are more difficult to 
pin down. Within the errors, two large sets of conditions 
yield line ratios consistent with the observed values. The 
first set is very similar to the arm conditions found above: 
about two thirds of the emission must come from relatively 
hot, fairly tenuous gas (Tki n ~ 100 K; rin 2 « 100 cm-' 3 ) 
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while the remainder is contributed by denser and cooler 
gas (T kin = 40±20 K; n Ha « 1000 cm" 3 ). The second set 
of possible solutions requires all emission to come from 
moderately dense (n H2 = 1000 — 3000 cm~ 3 ) gas. Half of 
the emission is then contributed by fairly cold (T kin sa 20 
K) and half by fairly warm (Tk; n ~ 150 K) gas. The 
present observations do not allow us to distinguish be- 
tween these two possibilities; both are equally likely. 

4.2. Beam-averaged molecular gas properties 

The chemical models presented by van Dishoeck & Black 
(1988) show a strong dependence of the N(G)/N(CO) col- 
umn density ratio (i.e. how much more carbon there is 
than the fraction tied up in carbon monoxide) on the total 
carbon and molecular hydrogen column densities. It thus 
provides a relation between the amounts of carbon monox- 
ide, neutral carbon and ionized carbon that produce the 
observed line intensities from CO, C°, and C + . In our 
analysis, we have assumed that the kinetic temperatures, 
H 2 densities and filling factors implied by the CO analysis 
equally apply to the [CI] and [CII] emission. In principle, 
one can then solve for column densities N(Cl)/dV and 
iV(CII)/dV". In practice, the column density of one com- 
ponent is usually well-determined, but that of the other 
is more or less degenerate. For this reason, we solved for 
identical velocity dispersions in the two gas components. 
Finally, we related total carbon (i.e. C° + C + + CO) col- 
umn densities to molecular hydrogen column densities by 
using an estimated [C]/[H] gas-phase abundance ratio. 

Although the analysis in terms of two gas components 
is rather superior to that assuming a single component, 
it is still not fully realistic. For instance, the assumption 
of identical beam-filling factors or identical velocity dis- 
persions for the various species ( 12 CO, 13 CO, C°, and 
C + ) is not a priori plausible. Fortunately, these assump- 
tions are useful but not critical in the determination of 
beam-averaged parameters. If, by way of example, we as- 
sume a smaller beam-filling factor, the model cloud in- 
tensity increases. This generally implies a higher model 
column-density which, however, is more strongly diluted 
by the beam. The beam-averaged column-density is mod- 
ified only by the degree of nonlincarity in the response of 
the model parameters to a change in filling factor, not by 
the magnitude of that change. 

Oxygen abundances and gradients for M 51 were taken 
from Vila-Costas & Edmunds 1992 and Zaritzky et al. 
1994. From the results published by Garnett et al. (1999) 
we estimate that at the relevant high metallicities [C]/[0] 
w 1, leading us to adopt a carbon abundance [C]/[H] = 
[0]/[H] = 3.0±0.5 x 10~ 3 . As a significant fraction of car- 
bon is tied up in dust particles and thus unavailable in the 
gas-phase, we have adopted a fractional correction factor 
S c = 0.27 (see for instance van Dishoeck & Black 1988), 
so that N H /N C = [2N(R 2 ) + N(BT)}/[N(CO) + N(CIT) 
+ iV(CI)] = 1250, uncertain by about a factor of two. 
In Tabled we present beam-averaged column densities for 



both CO and C (C° + C + ), and H 2 column densities de- 
rived under the assumptions just discussed, as well as the 
face-on mass densities. As the observed peak CO inten- 
sities are significantly below the model peak intensities, 
only a small fraction of the (large) beam surface area can 
be filled with emitting material. For position A (central 
region) we find a beam area filling factor of 70 ± 15 for 
either model, whereas we find a filling factor of about 20 
for position B (the spiral arm CO complex). 

In order to gauge the reliability of these results, we 
have explored the effect of the assumptions discussed 
above. If we do not assume identical velocity dispersions 
for the two gas components, but e.g. dispersions propor- 
tional to the kinetic temperature, derived H 2 columns and 
mass-densities change by amounts varying from 10% to 
20% in the three models described. If we do not assume 
equal beam filling factors for CO and C + (which greatly 
exceeds the contribution by C°, see below), but vary that 
of C + by e.g a factor of two, resulting H 2 columns and 
mass-densities change by amounts varying from 25% to 
30%. Differences in beam-filling factors cannot be very 
much larger, because then the CO and [CII] intensities, 
the [C]/[H] and [C]/[CO] abundances, and the H 2 column 
densities become mutually incompatible. 

4.3. Discussion 

As mentioned above, conditions in the center are fully 
described by either of two models that cannot be distin- 
guished as both are fully consistent with the observed line 
ratios and intensities. Observations with much higher res- 
olution using the Smithsonian Millimeter Array (SMA) 
recovered only a small fraction of the total J=3-2 12 CO 
flux. They showed, however, the existence of a compact 
(a few arcsec in size) molecular cloud centered on the nu- 
cleus (Matsushita et al. 2004) not resolved in our maps. 
For this cloud, the authors estimated high densities and 
temperatures that seem more or less consistent with those 
of our Model 2/Component 1. However, as they used only 
single-component LVG fitting such resemblance may be 
fortuitous. 

Our two models have in common that they explain the 
observed emission by requiring the existence of at least two 
gas components at different densities and temperatures. 
Nevertheless, they imply very different ISM conditions as 
these kinetic temperatures and densities vary significantly 
between models. The implied hydrogen column densities, 
mass surface densities and CO-to-H 2 conversion factor X 
also differ significantly, by a factor of three. That differ- 
ence is almost entirely caused by the fact that at lower 
kinetic temperatures a much larger ionized carbon (C + ) 
column density is required to produce the same observed 
[CII] line flux than at higher kinetic temperatures. 

Beam-averaged neutral carbon to carbon monoxide 
column density ratios range from N(C°)/N(CO) — 0.4- 
1.0, and are very comparable to those found in other 
galaxy centers (White et al. 1994; Israel et al. 1995; Stutzki 
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et al. 1997; Petitpas & Wilson 1998; Israel & Baas, 2001, 
2003). 

Our model results agree well with the parameters that 
we would have obtained by applying the PDR models pre- 
sented by Kaufman et al. (1999) to the observed intensities 
and ratios of the CO, [CI], [CII] lines, and the intensity 
of the far-infrared continuum (Rand et al. 1992). For the 
observed values, the PDR models predict the presence of 
molecular gas at temperatures Tkin ~ 100 K, and at gas 
densities 10 3 cm~ 3 illuminated by a fairly modest ambient 
radiation field with log G = 1.0 (G expressed in units 
of the Habing 1968 field: 1.6 x 10~ 3 erg s" 1 cm" 2 ). 

Our models also imply the central regions of M 51 to 
have a CO-to-H 2 conversion factor X = 0.50 ± 0.25 x 10 20 
H 2 mol cm~ 2 / K kms -1 , i.e. roughly a factor of four lower 
than the commonly used Galactic X value for the Solar 
Neighbourhood. However, such low values are not uncom- 
mon for galactic centers (Israel & Baas 2003 and references 
therein). For M 51, virtually identical results have also 
been obtained independently by Guelin et al. (1995) from 
an analysis of the 1.2/xm emission from cold dust, and by 
Nakai & Kuno (1995) when one considers their results for 
HII regions with one arcminute from the nucleus. Thus, 
we consider that the CO-to-H2 conversion factor X is un- 
usually well-established in the center of M 51; there can 
be little doubt that it is indeed much lower than the 'stan- 
dard' Milky Way value (see also Strong et al. 2004). We 
emphasize that low values for X mean that for any given 
CO intensity, beam-averaged molecular hydrogen column 
densities N( H2) and indeed H2 masses are much less than 
estimated from the 'standard' conversion. This has an- 
other important implication. All molecular species abun- 
dances derived by relating an independently determined 
column density to a molecular hydrogen column density 
obtained from the 'standard' Galactic conversion factor 
are in error and must be revised upwards by a potentially 
large factor. 

5. Conclusions 

1. Observation of the central arcminute of the well- 
studied spiral galaxy M 51 in various transitions of 
12 CO and 13 CO, and in [CI] shows that the molecular 
gas resides mostly in the bright inner spiral arms. 

2. The absence of a circumnuclear disk or torus in images 
of the center of M 51 whereas position- velocity maps 
exhibit a compact central feature exhibiting a rapid 
solid-body rotation implies that in other galaxies the 
presence of the latter cannot be taken as proof for the 
presence of the former. 

3. In the center of M 51, a dense component with 
n(H 2 ) ~ 10 3 cm -3 and T^n ~ 100 K is accompa- 
nied by either a less dense component with n(H2) ~ 
10 2 cm~ 3 at the same temperature, or a more dense 
component with n( H 2 ) ~ 3 x 10 3 cm -3 at a much lower 
temperature of Tkin ~ 20 K. Only a small fraction of 
all carbon is in CO. Total carbon column densities are 
about 7 times the CO column density. At least in the 



center of M 51, most of the carbon appears to be in 
the form of C + . The observed line ratios are consistent 
with standard PDR models. 

4. Emission from the CO peak in the northwestern arm 
likewise originates in (at least) two different gas com- 
ponents. One is modestly dense (n(H2) ~ 10 3 cm~ 3 ) 
and the other is relatively tenuous (n( H 2 ) « 10 2 cm~ 3 ) 
but more widespread. Both must be at elevated tem- 
peratures T Un = 100-150 K. 

5. The center of M 51 has a face-on gas mass density of 
40±20M Q pc~ 2 , and a relatively well-established CO- 
to-H 2 conversion factor X = 0.50 ± 0.25 x 10 20 H 2 mol 
cm _2 /Kkms _1 , i.e. a factor of four to five lower than 
commonly assumed standard Galactic values, but in 
line with similar determinations in other galaxy cen- 
ters. 
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